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Abstract 

Classical black holes and event horizons are highly non-local objects, defined in terms of the 
causal past of future null infinity. Alternative, (quasi)local definitions are often used in mathemat- 
ical, quantum, and numerical relativity. These include apparent, trapping, isolated, and dynamical 
horizons, all of which are closely associated to two-surfaces of zero outward null expansion. In this 
paper we show that three-surfaces which can be foliated with such two-surfaces are suitable bound- 
aries in both a quasilocal action and a phase space formulation of general relativity. The resulting 
formalism provides expressions for the quasilocal energy and angular momentum associated with 
the horizon. The values of the energy and angular momentum are in agreement with those derived 
from the isolated and dynamical horizon frameworks. 
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I. INTRODUCTION 



The basic ideas of black hole physics are well-known. These foundations, including the 
singularity, uniqueness, and area increase theorems pj, the laws of black hole mechanics 
0, and the identity between those laws and the laws of thermodynamics 0,0], have been 
understood for several decades. Perhaps the most remarkable conclusion to be garnered from 
these results is the realization that black holes, in spite of their apparently exotic origin, 
are not really so different from other objects studied by physics. Though classically they 
harbour singularities in spacetime and have interiors that cannot causally effect the rest of 
the universe, they still have physical properties such as mass, charge, angular momentum, 
temperature, and entropy. Nevertheless there is at least one fundamental way in which they 
are very different. With the exception of spacetimes that are known to be globally stationary, 
the exact location and size of a black hole cannot be determined by local measurements; 
instead the existence of a black hole event horizon is a property of the causal structure of 
the spacetime. Further, in traditional black hole physics, physical properties such as mass 
cannot be determined by local measurements, but instead are only defined with reference to 
asymptotic structures. This, of course, contrasts quite strongly with most physical objects 
which can be identified, measured, and understood without reference to the far future and 
spatial infinity. 

Given the non-local nature of event horizons, there has always been an interest in alterna- 
tive, local, characterizations of black holes. The original alternative black hole signature was 
the trapped surface [H[ and these figure prominently in the definition of apparent horizons 
PJ. More recently there have been several other (closely related) approaches. These include 
trapping 0, 0] , isolated |j| , dynamical [j| , and slowly evolving horizons . All of these 
will be discussed in more detail in Section Til Al but for now the important point to note is 
that they are all closely related to each other and defined by the convergence/divergence 
of particular sets of null geodesies. Specifically they all share one common property : the 
three-surface which they identify as a "horizon" can be foliated with spacelike two-surfaces 
of vanishing outward null expansion. These alternative notions of horizons are especially 
suited to studies of dynamic black holes. For an evolving black hole, the location of an 
event horizon is impossible to determine unless one first understands the causal structure 
of the entire spacetime. The alternative horizons do not face this problem as their location 
depends only on local geometry — the calculations involved may be difficult, but at least 
they are possible. 

Thus, alternative horizons are closely associated with the study of black hole evolutions 
and their interactions with their surroundings. The exchange of energy and angular momen- 
tum between a black hole and its environment is one of the most basic physical problems 
that might be investigated under these circumstances and this question has been addressed 
under each of the approaches. Various formulae have been proposed for energy and angular 
momentum associated with a horizon, as well as the flux of these quantities into and out of 
the hole. 

In this paper, we apply Hamiltonian methods to gain further insight into the energy and 
angular momentum of horizons. Specifically, we examine quasilocal action and Hamiltonian 
formulations of general relativity on manifolds which have horizons as boundaries. We begin 
in Section IIIII with a discussion of the action in the spirit of the well-known Brown- York 
quasilocal action formalism 11, ^| and its non-orthogonal generalizations (13 . 14j . We show 
that any three-surface (regardless of signature) which can be foliated with spacelike two- 
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surfaces of vanishing outward null expansion is a suitable boundary in a quasilocal action 
formulation for general relativity, provided that we fix the intrinsic geometry of those two- 
surfaces along with the connection on their normal bundles (a quantity closely associated 
with angular momentum). This boundary information only needs to be fixed in a weak form 
— up to diffeomorphisms which may "move it around" . 

If we impose a notion of time on the quasilocal region, then we can perform a Legen- 
dre transform of the action and so obtain a Hamiltonian for the spacetime. As in other 
derivations of a Hamiltonian, the bulk part of this functional is made up of constraints and 
so vanishes on solutions of the Einstein equations. Thus one can evaluate it on a slice of 
a solution by examining only the geometry of the boundary. Consequently, we can derive 
a quasilocal energy expression for these alternative horizons. As usual, the value of the 
Hamiltonian is sensitive to the precise form of the action. With our boundary conditions, 
there is a class of acceptable action functionals; the freedom in the action amounts to being 
allowed to add free functionals of the fixed data to the original form. As a consequence, the 
formalism fixes the energy only up to free functionals of the fixed data. 

In Section I1V1 we turn from Lagrangian to Hamiltonian arguments and consider the 
horizons as boundaries from the point of view of phase space, a symplectic structure, and 
Hamiltonian evolutions of a quasilocal region of spacetime. This investigation helps to 
further elucidate the allowed freedom in the definition of an energy and also gives us a 
concrete expression for the angular momentum associated with the horizon. The calculations 
follow the extended phase space formalism of 0] which is designed for situations where there 
are fluxes of quantities such as energy and angular momentum through the boundaries. This 
is in contrast to the standard Hamiltonian formalism which is designed for isolated systems 
where there can be no flux of these quantities (and so they are conserved). 

Applying this analysis, we derive expressions for the energy and the angular momen- 
tum associated with the horizon. The energy expression so derived matches the one found 
using the action/Legendre transform methods - essentially any diffeomorphism invariant 
functional of the intrinsic boundary two-geometry and the connection on its normal bundle. 
The new piece of information is the expression for the angular momentum. 

Finally, in Section we compare the Hamiltonian expressions for energy and angular 
momentum with those arising from the isolated and dynamical horizon calculations. On 
sections of the horizon which are null, we find that the angular momentum is identical to 
that obtained in the isolated horizon framework. Additionally, if we restrict our generators 
of rotation to divergence free vector fields, the associated angular momenta on spacelike 
sections of the horizon agree with the dynamical horizon values. Next, we turn to the 
energy. We can greatly restrict the freedom in the energy by requiring that it satisfy a 
basic property: the value of the energy must transform appropriately with rescalings of 
the evolution vector field. We find that in this case the natural evolution vector field is 
the null vector, £, for which the outward expansion vanishes. This is true even though it 
is not, in general, tangent to the horizon. Furthermore, the energy associated to a given 
evolution is then determined up to an overall constant rescaling. This constant can be fixed 
by appealing to the dynamical horizon framework. With this choice, the energy derived 
using Hamiltonian methods will satisfy the dynamical horizon flux law. 
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II. LOCAL HORIZONS - A REVIEW 



A. Definitions 

There are several (quasi-) local definitions of a horizon. In this section, we summarize 
the various definitions, highlighting both their similarities and differences, and emphasizing 
those features relevant for the action and Hamiltonian calculations that will follow. 

Spacelike two-surfaces and their null normal vectors feature prominently in each of the 
alternative horizon definitions, so we begin by discussing them. For any spacelike two-surface 
S, we can find families of null vectors that are normal to the surfaces. Up to normalization 
and orientation, there are just two such null directions. If the spacetime is time-oriented 
we can single out the forward-in-time pointing directions and so, up to rescaling, define two 
future-pointing null normals £ a and n a . Then, if q ab is the two-metric on S, the expansion 
associated with these vector fields is given by: 



These are expansions in the sense that if one Lie-drags S in the direction £ a , and e is the 
area two-form on S, then 



Thus, 9<£\ describes how the area-form on S expands/contracts along £. Similar arguments 
apply to 9 {n) . 

The two-surface 5* is said to be trapped if it is closed and the expansions of both £ and n 
are negative at every point on S. A marginally trapped surface is a slight generalization of a 
trapped surface; in this case, the expansions of the null vectors are required to be less than or 
equal to zero. Trapped and marginally trapped surfaces can be contrasted to surfaces in flat 
spacetime where one would expect the outward expansion of convex surfaces to be positive 
everywhere, and the inward expansion negative. Intuitively, the negative expansion of both 
null normals can be attributed to strong gravitational fields in the interior of S. This picture 
is strengthened by the singularity theorems which state that (provided energy conditions are 
satisfied), any spacetime containing a trapped surface must also contain a singularity [Hll6j|. 
Furthermore, it can be shown that any trapped surface in an asymptotically flat spacetime 
must be contained within an event horizon. Even in the absence of an appropriate asymptotic 
structure, a trapped surface may always be extended into a spacelike three-surface that hides 
its inside from its outside [T3|. Thus, it seems natural that the various notions of a horizon 
discussed below all make use of trapped or marginally trapped surfaces. 

We begin with a brief discussion of apparent horizons Let £ be an asymptotically 
flat Cauchy surface. From an intuitive point of view, such a £ is an "instant" in an asymp- 
totically flat spacetime. Now, a point p G £ is said to be trapped if it lies on a trapped 
surface that in turn is contained entirely in £; the total trapped region T is the union of 
all the trapped points in E. An apparent horizon is then defined as the boundary A of the 
total trapped region T. With certain assumptions about smoothness ^ij], it can be shown, 
as one would expect, that 9m = on A . Further, it is true that T is itself contained within 
an event horizon. In the case of stationary black hole solutions, such as Schwarzschild, T 
completely fills the region of £ inside the black hole and the apparent horizon coincides with 
the event horizon. 



00 = q ab V a £ b and 9 {n) = q ab V a n b 



(2.1) 



(2.2) 
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While it is certainly possible to identify an apparent horizon without reference to the far 
future, it is not a trivial task to first locate all the trapped surfaces, find their union and then 
obtain the boundary of that region. Further, one requires knowledge of the entire "instant" 
before one can locate the apparent horizon. Finally, A won't always satisfy the smoothness 
assumptions mentioned in the previous paragraph (see [3] for a good discussion). All of 
these factors can make working with A inconvenient (at the very least). 

More recent quasilocal horizon definitions attempt to extract the idea of apparent horizons 
as the boundaries of regions of especially strong gravitational fields, while dropping some of 
their more unwieldy aspects. Specifically they circumvent a consideration of entire trapped 
regions and focus directly on two-surfaces S of vanishing outward null expansion 9m = 0. 
This has the advantage of making the definitions quasilocal. Furthermore, as is the case for 
apparent horizons, in asymptotically flat spacetimes such surfaces S are always contained 
within event horizons. In common usage, for example in numerical relativity |20j, the term 
"apparent horizon" is frequently used to refer to the outermost two-surface for which 9m = 
as it is simply not practical to find the boundary of the total trapped region. However, one 
should keep in mind that while in some circumstances the definitions are equivalent, these 
surfaces do not always coincide with A. 

Notions such as trapping, isolated, and dynamical horizons seek to formalize these vernac- 
ular ideas of an apparent horizon and at the same time remove their somewhat problematic 
dependence on £ slicing. 1 Rather than considering two-dimensional surfaces with vanishing 
outward expansion in some three-surface S, they instead define three-dimensional horizons 
within the full four-dimensional spacetime. These horizons are required to be foliated by 
closed spacelike two-surfaces S v on which 9m = 0. Thus, while slicings of the spacetime 
may be found that are compatible with these 9m = surfaces (in which case they will often 
be apparent horizons), the horizon definitions are independent of any such foliation of the 
spacetime. 2 

We will begin by considering isolated horizons and then move on to the very closely 
related trapping and dynamical horizons. For detailed discussions of the merits of each 
approach the reader is referred to the literature. Here we will briefly review their relevant 
properties and demonstrate that each falls within the scope of the upcoming calculations. 

Isolated horizons (if are intended to characterize black holes that are not interacting 
with their surroundings, even though those surroundings may themselves be dynamic. An 
isolated horizon is a null three-surface whose null normal has zero expansion. Every foliation 
of an isolated horizon will have the same outward null normal, whence the foliation into two- 
surfaces is not important. Making use of a mild energy condition, it is straightforward to 
demonstrate that isolated horizons don't change in either shape or size and no energy flows 
across them. With an additional restriction on the permitted class of null normals £ a , it 
can be shown that isolated horizons obey the zeroth law of black hole mechanics, namely 
that the surface gravity is constant over the horizon. Furthermore, a phase space version 
of the first law has been shown to hold. Finally, loop quantum gravity calculations have 
demonstrated that they have an entropy proportional to their surface area 2^ . 



1 For example, it is well-known that even Schwarzschild spacetime may be foliated so that there are no 
apparent horizons in any slice |2l| . 

2 Note however that the slicing dependence discussed in the last footnote is then replaced by a non-rigidity 
of the horizons. In general the two-surfaces making up these horizons may be smoothly deformed while 
keeping their 9^ = property. See for a recent discussion of this non-uniqueness. 
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Trapping horizons p are defined as three-surfaces (of indeterminate signature) that can 
be foliated with 9m = two-surfaces S. By themselves these conditions are not enough to 
distinguish between outer and inner black hole horizons or for that matter between black 
holes and white holes. Additional conditions are added to define future outer trapping 
horizons (FOTHs). A FOTH is a trapping horizon that also satisfies the property that 
9(n) < and C n 9^ < 0. Thus, ingoing null geodesies are converging, while outgoing 
null geodesies converge inside the horizon, diverge outside and have zero expansion on the 
horizon. Given a three-surface E that intersects a FOTH along a foliating two-surface S, 
these conditions are intended to ensure that the region inside of S is trapped, while the region 
outside will not be (thus S should be apparent horizon). In a spacetime respecting the null 
energy condition, the above conditions are sufficient to guarantee that dynamic trapping 
horizons are spacelike and expanding, and so individually obey a form of the second-law 
of black hole mechanics. In the absence of (gravitational wave) shearing or matter falling 
across the horizon, they are null and non-expanding. Finally for energy condition violating 
matter, horizons can shrink in size, in accord with Hawking radiation. 

A slowly evolving horizon 0] is an "almost" isolated FOTH and so defined in the same 
way but with an extra condition that characterizes "almost". Intuitively, these horizons can 
be thought of as FOTHs which are only weakly interacting with their environment. In this 
regime, a truly dynamical version of the first law for black holes has been formulated. This 
is also the regime of interest for studying the approach to equilibrium after a gravitational 
collapse or black hole merger. 

A dynamical horizon (DH) is also closely related to a FOTH. It is a spacelike three- 
surface that can be foliated with spacelike two-surfaces for which < and 9^ = 0. Here, 
the requirement that the surface be spacelike has replaced the C n 9(e) < condition of the 
FOTH. This ensures that the horizon can be identified without the need to take derivatives 
out of the three-surface. In cases where the horizon is expanding (which are largely the cases 
of interest) the Einstein equations along with energy conditions ensure that the definitions of 
dynamical horizons and FOTHs are equivalent. It has been shown that dynamical horizons 
satisfy a variety of conservation laws that relate the change in horizon properties to fluxes 
of stress-energy and matter across the horizon. The connection between those laws and the 
calculations of this paper will be discussed further in Section IV! 

Finally, there has recently been some discussion of non-FOTH/DH extensions of 
FOTHs/DHs HI H @. Referred to as marginally trapped tubes (MTTs), these are 
three-surfaces that satisfy 9(e) = and 9^ n ) < but are not not necessarily everywhere null 
or spacelike. In the terminology of trapping horizons, they are future trapping horizons 
with the slight generalization that £ n 9(e) may take any value. It has been shown that a 
single MTT can include not only regions that are isolated and dynamical but also regions 
with timelike signatures (where £ n 9(e) < 0). Such regions appear to be associated with 
extreme conditions where new horizons form outside of old ones. Though these so-called 
timelike membranes probably cannot be viewed as horizons (they are traversible in both 
directions), they can evolve directly into/from DHs and so need to be kept in mind during 
the following discussions as possible quasilocal boundaries. That said, MTTs will usually 
be excluded from consideration in the following as they introduce complications such as 
non-trivial intersections with any foliation Ej of the spacetime (see [24| for many examples 
of these complications). 
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B. Horizon Geometry 



In the preceeding section, we have discussed many different local horizon formulations. 
Here, we would like to examine the geometrical properties of the horizons, focussing in 
particular on those which are common to all the formulations. Initially, we will restrict 
our attention to two-surfaces S of zero expansion which foliate a three-surface H. The 
intrinsic geometry of the horizon is completely determined by the two-metric q ab and the 
corresponding area two-form e introduced previously. The extrinsic geometry of the surface 
is partially characterized by the extrinsic curvatures (k ab ) of S along the null directions, 

2%)?a& + QacQbdrfj!) = k {e)ab := q a c q b d VJ d and 

^0(n)Qab + qacQbd^) = k {n)ab : = q a c q b d V ' c n d . (2.3) 



Here, 9^ and 9( n ) are the expansions of the null normals, introduced previously, while afy 
and a?K are the shears along £ and n respectively. Since £ and n are null normals to the 
surface S they are, by definition, twist free. The expansion of £ is required to vanish on all 
the horizons discussed above. The other quantities will generically be non-zero except on an 
isolated horizon when the shear afy must also vanish as there can be no flux of gravitational 
energy across the horizon. 

The rest of the extrinsic geometry is described by the connection on the normal bundle 
to S, 

uJ a = -n c V a £ c . (2.4) 

This quantity is closely associated in many of the formalisms with physical characteristics 
of the horizons. On an isolated horizon, the surface gravity is obtained by contracting with 
£, whence Kin = £ a uj a = —£ a n c V a £ c . Additionally, the two-surface components of u), namely 

Cu a = q a b u b = -q a b n c \7 b £ c , (2.5) 

are often associated with quasilocal angular momentum |H |j| [ll[ HfJ |2?J . Specifically, on 
an axisymmetric isolated or slowly evolving horizon with a symmetry generated by a , the 
angular momentum is given by 

■W] = -r?7 / e^a- (2-6) 

8ttG Js 

The definition of angular momentum on a dynamical horizon is slightly different. Since 
dynamical horizons are by definition spacelike, we can introduce the timelike unit normal 
to the horizon f a and the unit spacelike normal to the leaves of the horizon, f a . Then, the 
angular momentum is constructed from these vector fields as 

o-£^b\ 

tip 

S 



JdM = s^g ' tPrVah. (2-7) 



Then, it is straightforward to integrate the diffeomorphism constraint over the horizon to 
obtain a flux law for angular momentum . Specifically, given a region of the horizon AH 
bounded by two surfaces 52 and Si, the angular momentum satisfies 

jm - Ji[<t>\ = fj { T ^y b + i^g p ? c ^} ■ ( 2 - 8 ) 
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In the above e is the natural volume form on the horizon. The two terms in ()2.8|) can 
then be interpreted as the flux of matter and gravitation angular momentum through the 
horizon. 

III. THE ACTION 

In this section we construct an action formulation for general relativity over a quasilocal 
region of spacetime that includes a trapping horizon H as one of its boundaries. To do this, 
we first introduce the spacetime M of interest, and then we specify appropriate boundary 
conditions on each boundary of the spacetime. Those on the horizon will be specified 
in terms of the geometrical properties discussed in Section III Bl while those on the other 
boundaries are standard and well-known. Combining these with an action whose bulk term 
is the usual Einstein-Hilbert term, we show that the first variation of the action vanishes 
in the required way if the Einstein equations hold in the bulk. Thus, an action formulation 
exists for spacetimes which have a horizon as a boundary. Finally, by performing a Legendre 
transform, we derive a corresponding quasilocal Hamiltonian from this action. In the usual 
way, the form of this functional then provides a definition of the energy associated with the 
region and its assorted boundaries. 

A. The spacetime and its boundaries 

Let M be a finite region of spacetime with metric g a b, metric-compatible derivative V a , 
and four- volume form tahcd (also written as e). We require M to be time orientable and 
diffeomorphic to [0, 1] x [0, 1] x S 2 , where the four boundaries are: 

1. a spacelike past boundary Si which represents the initial configuration, 

2. a spacelike future boundary S 2 which represents the final configuration, 

3. a timelike outer boundary B, and 

4. an inner horizon boundary H which is foliated by two-spheres 3 S v of zero outward 
null expansion {6m = 0). 

The two-dimensional intersections of these boundaries are B± — Si n B, £>2 = S 2 n B, 
Si = Si H H, and S2 = S2 D H. Figure ^ provides a three-dimensional representation of the 
manifold M and these assorted boundaries. 

In describing the precise boundary conditions on each of the surfaces and also in the 
ensuing calculations, it will often be important to distinguish between four-dimensional 
quantities and those defined exclusively in the various two- and three-dimensional surfaces. 
To this end, we introduce some notational conventions. Tensor fields which are elements of 
the tensor bundles of the full spacetime will be labelled with lower case early alphabet Latin 



3 Strictly speaking the only property of S 2 that we will need in the following calculations is that it is closed. 
That said it can be shown [|J that in most physically interesting situations, the Einstein equations along 
with energy conditions actually restrict the S v to being topologically S 2 . Thus, for definiteness, we make 
this assumption. 
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B 2 




FIG. 1: A three-dimensional representation of M along with its boundaries. 

indices (a, b, c, etc.). Tensors in the bundles associated with three-dimensional hypersurfaces 
will be labelled with lower case mid-alphabet Latin indices (i, j, k, etc.) and tensors in the 
bundles associated with two-dimensional surfaces will be labelled with upper case Latin 
indices (A, B, C, etc.). Additionally, we will need to map tensors between the tensor bundles 
associated with the different manifolds. We will write all pull-back/push- forward operators 
as p with indices added to indicate which spaces it maps between. Thus, the four-dimensional 
metric g a b pulled back into a two-dimensional surface would be written PAPsdab, while a 
three-dimensional vector v l would push forward to p^v 1 in the full spacetime. Of course, 
as indicated by the indices, only fully covariant tensors may be pulled back and only fully 
contravariant tensors may be pushed forward in a metric independent way. If we wish to 
project tensors of a different rank we must make use of the metric, rendering the results 
metric dependent. 

With the preliminaries out of the way, we now describe in more detail the various bound- 
aries of M. Figure El provides a two-dimensional representation of the spacetime M and its 
boundaries. Surface normals and other important vector fields are also shown. 

Initial and final boundaries : Si, £2 

Denote the future directed unit normal vectors to the spacelike boundaries £^2 by u a . 
Then, the three-metric on the surfaces E is given by hij which satisfies h ab := p^p^ = 
gab _|_ ^a^& me tric compatible derivative operator is Di. We obtain the three-volume 
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Si 



FIG. 2: A two-dimensional representation of the spacetime containing M (the lightly shaded 
region), showing the past and future three-boundaries Si and £2 along with the inner and outer 
three-boundaries H and B. The intersections of these surfaces Si 2 and B± t 2 are also shown along 
with the orientations of the assorted normal vector fields for the case where H is spacelike and 
expanding. The darkly shaded region is the interior of H - the "hole" . Each point of this schematic 
is a two-sphere in the real M. 

form by contracting the four-volume form with the normal u a . A Thus, (E e a fe c = e a b c dU d from 
which it follows that e = u A e. Finally, the extrinsic curvature of E in M is given by 
Kfy = h ac h bd V c u d . 5 

Outer boundaries : B, Bi,&2 

In keeping with the standard quasilocal action analyses, we take the outer boundary B 
to be timelike with n a as its outward pointing spacelike unit normal. Then, the three- 
metric on B is 7^, where 7°* = p i a Pj b, ~f^ — 9 ab — n a n h and it defines the metric compatible 
derivative A$. The extrinsic curvature K?^ = / -f aC/ -f bd 'V c nd and the volume form is taken to 

be {B e abc = -e abcd n d , so that e = n A (S e. 6 

The outer boundary B intersects Ei and E 2 on the spacelike two-surfaces B± and &2- 
We do not require that this intersection be orthogonal (i.e. there is no requirement that 
u-n = 0). Thus, we introduce two additional normals: n a , the outward unit spacelike normal 



4 The volume form is induced on every submanifold by Stokes theorem. It can be found by contracting 
Cohcd with the outward normal on the last index. 

5 For those familiar with the work of Brown and York , note that they define extrinsic curvature with 
the opposite sign. 

6 Note that this is the opposite of the volume form induced by Stokes theorem. Consequently we will 
introduce a minus sign when integrating an exact derivative out to the boundary B. 
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to Bip in £1,2, and u a , the future directed timelike unit normal to /3x,2 in B (as illustrated 
in figure EJ) • Defining 77 = —u ■ n, the unit normals are related by 

u a = y/l + r] 2 u a + r]h a and n a = 77 u a + y^T+rf h a . (3.1) 

Clearly, when 77 = 0, u a = u a and n a = n a and so the H/B intersection is orthogonal. 

The two-metric on /3 12 is denoted jab and its inverse pushes forward to 7 afc = h ab —h a n b = 
7 ab + u a u b . The derivative operator compatible with this metric is d^. The volume form e 
on #1,2, is given bye = uAnAe = uAnAe. 

Inner boundaries: H , Si, S2 

A discussion of the inner boundary H is more involved than that of the outer boundary. 
Since this boundary is required to be a trapping horizon: 

1. H has additional structure in the form of a foliation into the 9m = spacelike two- 
surfaces and 

2. there is no restriction on the signature of H and so it can be spacelike, null, timelike, 
or conceivably some mixture of the three. 

We label the foliation two-surfaces S v , where v is a foliation parameter. Furthermore, we 
assume that Si = Ex Pi H and S2 = E2 H H are the initial and final leaves of that foliation. 
That is, Si = S Vl and 5* 2 = S V2 and 7J1 < v < v<i. The foliation parameter naturally induces 
the one-form dv G T*H normal to the horizon 2-surfaces. In addition, we introduce a vector 
field tangent to the horizon and satisfying C(q/q v )V = 1. It is not unique, but will be 
useful in describing how the cross-sections S v evolve. 

Since the horizon is foliated, we can introduce normals to the two-surfaces. It is most 
convenient to work with the null normals l a and n a to the surfaces S v . This will allow us to 
treat timelike, spacelike and null horizons H together and also make it easier to impose the 
condition that 9m = 0. As shown in figure El we let l a be outward (towards B) and future 
pointing while n a will be inward (away from B) and also future pointing. Furthermore, we 
require the standard normalization £-n = — 1. This, however, still leaves one rescaling degree 
of freedom in the null vectors which will be important during the calculations. Irrespective 
of how this is chosen, the two-metric (Jab satisfies q ab = p A a p^q AB = g ab + £ a n b + n a £ b . As on 
the other two-surfaces, the metric compatible derivative will be denoted 4a- Additionally, 
in Section Til Bl we have introduced the extrinsic curvatures and kfS of the null vectors, 
as well as the one form Co a which contains the angular momentum information. 

On the initial and final cross sections of the horizon, we can relate the null normals I and 
n to those inherited from Ex and E 2 . Specifically, we have the timelike normal u a inherited 
from the Ex/E 2 boundaries and the outward-pointing spacelike normal h a to Si/S 2 in Ex/E 2 
as shown in figure 121 Then 

t = ((u a + n a ) and n a = ^-{u a - h a ) , (3.2) 

for some function £. Furthermore, we can also write the two-metric as q ab = g ab +u a u b —h a h b . 

We now turn to characterizing the three-dimensional horizon H. We begin by introducing 
a normal to H in M, which we denote r a , and a vector field V a that is both tangent to H (so 
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that V a r a = 0) and perpendicular to the S v . Both are assumed to be future pointing (in the 
sense that V • n and r • n are strictly negative). Since the horizon is not of fixed signature, 
there is no natural normalization for these structures, however, it will often be convenient 
to normalize both relative to the foliation label v. Thus, we require 

d 

C v v = l ^ — = V + V where V G TS V . (3.3) 

ov 

and furthermore, 

r . r = -V-V. (3.4) 

For those familiar with dynamical horizons when H is spacelike, V a is spacelike and 
proportional to the unit vector f a of those papers. Similarly r a becomes timelike and is 
proportional to the unit f a of those papers. By contrast, when H is null, the fields are 
degenerate and are both parallel to £. 

For a given two-surface S v , we have introduced two sets of (future-pointing) normal vector 
fields: (£ a , n a ) and (V a , r a ). Clearly these must be related. To this end, we introduce a scalar 
field C on H chosen so that £ a + Cn a is normal to H, and hence proportional to r a . It follows 
that l a — Cn a will then lie in H and be proportional to V a . Therefore, we can write 

V° = a{t - Cn a ) and r a = a(£ a + Cn a ) , (3.5) 

for some a > on H. The sign of C naturally relates to the signature of the induced metric 
Qij = Pi a Pj b 9ab on H; specifically, if C < 0, C = 0, or C > 0, H is respectively timelike, null, 
or spacelike. Finally, the orientation of the volume form on the horizon is chosen so that 

£ = rAduA£ = -lAnAe. (3.6) 

Note that we have said very little about the geometry of H itself. The reason for this is 
that the ensuing calculations will largely focus on two-surface rather than three-dimensional 
quantities. This is the most convenient way to proceed given that H does not have a fixed 
signature and so may change between spacelike, null, and timelike signatures. Similarly, in 
order to avoid having to deal with the timelike, null, and spacelike quantities seperately we 
choose not to work with normalized quantities such as f a and f a which are ill-defined if H 
becomes null. 



B. Boundary conditions 

In order to obtain a well-defined action principle, it is necessary to specify boundary 
conditions for all of the boundary manifolds presented in the previous section. First, on the 
initial and final boundaries Si and S 2 , as well as the outer boundary B, we fix the three 
metric. Similarly on the two-dimensional intersections £>, we fix the two-metric. Thus, we 
require that 

Shij = , 5^ij = and S'Jab = . (3.7) 

These are by no means the only permissible boundary conditions, but they are by far the 
most common. See |l2| for a discussion of alternative boundary conditions. 

The situation on H is more complicated and consequently more interesting. We start with 
a defining condition, 9^ = 0, and so expect it to play an important role in any calculations 
- otherwise the resulting action formulation will in no way be tied to horizons. With this 
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principle in mind we reject the idea of fixing the full three-metric qij as this would give a 
well-defined action principle regardless of the value of 6^ . Geometrically, it would also seem 
to be overly profligate compared to what is done on the other boundaries. For example, for 
Ei, the geometry is fully specified by its intrinsic metric, hij, and extrinsic curvature, K?Ly 
Each of these has six components, so for Ei with standard boundary conditions, six out of 
a possible twelve geometric components are fixed. Now, if we fix q^ on H as well as 6^ = 
we would essentially have fixed seven of those components. By contrast, the path that we 
follow in this paper will fix 6m = along with just five other components of the surface 
geometry on H, and each of those conditions will play a crucial role in the calculations. 

H is also different from the other boundaries in that it comes with a foliation. This 
foliation is not just an optional structure imposed on H, but is crucial to its very definition. 
Without the foliation, the key dm — condition is not well defined and, furthermore, for a 
given (spacelike) H this foliation is unique j25[. Given this intricate association of H and 
its foliation, it is perhaps not surprising that it also features prominently in our boundary 
conditions on H. We fix: 

1. e (£) = o, 

2. the spacelike two- metric q~AB on each S v (which also fixes e), 

3. u a = —PA an b^a^ b (the connection on the normal bundle to each S v ), and 

4. ^ = Pi a £ a , the pull-back of i a to T*H. 

Then in the action formulation, the first variation of each of these quantities is zero. The 
foliation surfaces S v as well as their coordinate label v are also kept fixed under these 
variations. 

The first condition is simply one of the defining conditions for a trapping or dynamical 
horizon. Note however, that other defining conditions such as 6i n ) < 0, C n 9m < 0, or H 
being spacelike, are not imposed. Thus, these boundary conditions apply equally well to 
dynamical horizons or any of the varieties of trapping horizons, including those usually taken 
to represent white holes or inner black hole horizons 0. It also includes three-dimensional 
surfaces traced by smoothly evolving apparent horizons. 7 

The second condition fixes the intrinsic geometry of each of the S v two-surfaces and so 
is somewhat similar to the usual fixing of the boundary metric. The difference, of course, is 
that here only the two-geometries are fixed. The rest of the three-geometry which gives the 
details how the two-surfaces evolve into each other, including the signature of H, is left free. 
Thus, while we know quantities such as the curvature scalar and area of the two-surfaces, 
the boundary conditions won't fix such things as the shear of H, or indeed any other locally 
defined rate of change. 

The third condition fixes another part of the extrinsic geometry of S v by specifying the 
connection on the normal bundle. Equivalently, it may be thought of as fixing the angular 
momentum of the horizon. Combined with the previous conditions, we see that the boundary 



7 However, we do not allow for situations where C n 9(e\ changes sign and so C diverges, such as the space- 
times considered in Q. Those stations can be dealt with fairly easily by sectioning the horizon into 
spacelikc/timclike sections over which C is well-defined, but for simplicity we will not consider these cases 
here. 
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conditions on H are mixed in the sense that they fix parts of both the intrinsic and extrinsic 
(two) geometry of the S v . This contrasts with the conditions on B an d £i ? which exclusively 
fixed the intrinsic (three) geometry. It is interesting to note that Epp [27[ presented counting 
arguments which suggest that knowledge of the intrinsic geometry of a two-surface plus the 
curvature on its normal bundle should be sufficient to uniquely determine its embedding in 
a four- dimensional space of constant curvature. Thus, in this sense, the information that 
we have specified on S v can be thought of as a complete set of data that fully specifies its 
geometry. 

Finally, we come to the fourth term. At first glance it may appear that we are fixing 
a seventh "component" of the geometry (given that 9^ = is one, (jab gives three, and 
Co a gives two). However, that is not the case. Recall that when defining l a we left its 
normalization (and so that of n a ) unspecified. This condition removes that degree of freedom. 
Hence, it is essentially equivalent to a requirement like the one that specifies that extrinsic 
curvatures should be defined using unit-normal vectors to surfaces. Just as such a definition 
is necessary so that a discussion of extrinsic curvature is meaningful, here we need a fix to 
give Co a a geometric meaning — otherwise it would change with rescalings of £ a . Keep in 
mind that we only fix the pullback of £ to H, not i a itself (which would have geometric 
significance) . 

C. The quasilocal action and its variation 

Having defined M and established this notation, we turn to the action itself 
for the region M and its boundaries is given by 

Ih = TTTr l en + lTr 

16ttG J m 8ttG i E2 _ Sl 

-8^/j B ^ v - + 8^/ B :_ s f" (,,) 

where 1Z is the Ricci scalar over M, k\> is given by 

Kv = V a uo a = -V a n b Vj\ 

and V a is defined in ()3.3|) . Equivalently, Ky + V a Coa = (~§^Y u u an d so we see that this 
action explicitly depends on both the foliation and the associated vector field (J^) (though 
in the upcoming calculations, we will see that this is essentially a gauge dependence only as 
the variation of Iff will vanish whatever the value of these coordinate quantities). 

The bulk, £, B, and B boundary terms are standard and found in all discussions of the 
gravitational action. The integral of the Ricci scalar over M is the Einstein-Hilbert bulk 
term while the integrals over Si, £2, and B are the usual extrinsic curvature terms which 
must be included when M has timelike/spacelike boundaries (see, for example, 0]). The 
fourth term may be a little less familiar, but is also well known |2§j. This "corner" term is 
necessary when the £ boundaries do not orthogonally intersect B\ that is u a n a 7^ 0. The 
final two terms arise at the horizon boundary and so are new, though we immediately note 
that the last is a "corner" term between H and £1,2- 



The action 



(3.8) 



(3.9) 
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Evaluating the first variation of this metric is a non-trivial calculation. Most readers will 
not be interested in the details and so they are relegated to Appendix [X] Here, we just 
quote the final result: 



SI ° = T^g jf^ ~ itc IS<^ ~ KJ, " } Sh « (3 ' 10) 

[dvAe(S AB 5q AB ) + -^- [se(l-ln() 

J H 87TG Js 2 -S 1 

+ 8^G J H i* A {( 6 & 9 V> ~ + VA 6 ( dv ^^a)-5 (20 {e) n Ae)}. 

Here, the quantity s AB is given by 

s AB = (k T + K V )q AB -k? B , (3.11) 

while ^_ and n respectively denote the pull-backs of £ a and n a into H . 

As for the action itself, the boundary terms on Si 1 2 and B\ 1 2 are completely standard 
and it is clear that fixing Shij = and Sjij = will ensure that these vanish. Similarly, it is 
not hard to see that our boundary conditions from Section Qll Bl are sufficient to ensure the 
terms on H and Si 2 vanish. Thus, 51 h = if and only if Einstein's equations hold in the 
bulk and so this action is well-formulated. 

Finally we note that, as in most Lagrangian analyses, the action that we have considered 
is not the only action compatible with the given boundary conditions. In fact, any functional 
of the form 

l' H — Ih + (functional of the fixed data) , (3-12) 

will work equally well, since the variation of a functional of the fixed data is zero. The 
implications of this freedom will be discussed in more detail in Section 1111 El 

D. Weak boundary conditions 

In the previous section, we have seen that the action principle is well-formulated for a 
manifold with a trapping horizon boundary. In the calculation, the boundary conditions have 
been imposed in the standard way; namely by requiring that variations of fixed quantities 
vanish. However for boundary conditions such as those on B, it has recently been shown 
[Til 15| that a weaker form of the boundary conditions can be used instead. These weak 



+ 



1 
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boundary conditions only fix the data up to diffeomorphisms. The motivation for this is the 
principle that physical theories should be generally covariant. In the presence of boundaries, 
we take this to mean that while we may (partially) specify the geometry of those surfaces, 
we should not tie that geometry to specific points of Si j2 , B± t 2, <Si,2, B or H — just as 
in general relativity spacetime geometry isn't tied to points of an underlying manifold (or 
equivalently a coordinate system). 

We now show that the boundary conditions of Section IIIIBI can also be applied in this 
weak form — up to diffeomorphisms which preserve the underlying manifold structure of 
the boundaries. Therefore, we require that 

Shij = C-vhij , 5%j = Czjij and S^ab = £>ylAB (3.13) 
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where Y l e TE 1>2 , Z* e TB and y A G TB ly2 . Since the boundaries B± :2 must be preserved 
(and so mapped into themselves), we impose the further restrictions that 

Y a = Z a = y a on B ly2 , (3.14) 

where each of these is the push-forward of the appropriate lower-dimensional vector field. 

Similarly, at the horizon, we allow diffeomorphisms which preserve the geometric struc- 
ture, which in this case includes the foliation S v . Thus, such maps should send all points 
in a given leaf S v into a second leaf S^xH f° r some function 4> x (v) '■ [^1,^2] - ► [ v i, v 2\- The 
most general vector fields generating a one-parameter family of diffeomorphisms of this type 
take the form 

X i = x V i + pjx A where x A e TS V , (3.15) 
where x Q = x Q (v) is a function of the foliation parameter v, and 

x = and Y i = on S 1)2 . (3.16) 
Then, the weak form of the boundary conditions fixes : 

Sq AB = CxQab , Su A = C x u A , 8 gj = Cx£_i and 6» M = 59 {i) = . (3.17) 

The action principle will remain well-defined if we also enforce the diffeomorphism con- 
straint on each of the boundary three-surfaces. This is not a big assumption as these are 
naturally induced if the Einstein equations hold in the bulk. Given a normal Xa and a 
tangent vector X a to a surface, the (vacuum) diffeomorphism constraint on that surface is 
X a G a b Xb = 0. Then on Si 2 and B this constraint says that 

A W - ^uh ij ) = and A< (Kg - K nl ij ) = , (3.18) 

where we respectively take u and n as the normal vectors. 

The situation on if is a little more complicated since its uncertain signature means that 
we cannot utilize unit normals or assume the existence of a (unique) inverse metric. The 
most useful form for our purposes is found by taking r a (defined in equation ([3.5)1 ) as the 
normal and evaluating 

X a G a \ = 0, (3.19) 

where we continue to assume that T a t, = 0. For the current calculation we further specialize 
this by integrating over S v . Then, as noted in [Io| : 
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I [K V C x e + i A C v [eCu A ]} = ! e {k AB C x q AB + 2aCC x 9 [n) } . 



(3.20) 



A full derivation of this result will be found in |29j , but for now a few lines of algebra puts 
it into a more useful form for the current calculation: 

Je{x A u A ) = ^JjdvAe(^C x q AB ^j - e A (C x £_)6 {n) - V A C x (dv A ec^)j . 

(3.21) 



8ttG 
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With the weaker boundary conditions, the variation of the action (|3.10|) becomes 



f ISO (A « " ^ £ + 8^? Stah " 1( " ) 
l - / *, A i rf„) + ^ / c_ f i ( 1 - in 
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+ A X ^ (n) " V A C x (dv A i£ A ) - 5 (2 n A e) } . 

The diffeomorphism constraints of £1,2 and B allow us to rewrite the boundary terms on 
those surfaces as exact derivatives. These can then be integrated to give contributions on 
B± } 2 and S\^- Similarly, the horizon constraint (|3.21|) can be used to replace the terms on 
H with a term on Si 2 - Therefore, we arrive at 

5I h = TTTT^ f^G ab 6g ab + -^— [ {(£ y e) sinh" 1 ^) - ey a (n b V a u b + u b V a n b ) } 

+^ / {(jC^)(l-\nC) + ex a (n b V a u b + n b V a u b )} (3.23) 

Finally, we can make use of the relationship between 77, n, u and n, u given in ()3.1j) 
to show that the terms on £> 12 evaluate to an exact derivative which then integrates to 
zero. Similarly the terms on Sip are exact as can easily be seen by using the relationship 
between (, n, u and £, n given in (|3.2|) . The only remaining term is then the bulk term which 
vanishes for solutions to Einstein's equations. Thus, we have shown that the action principle 
remains well formulated when we enforce weak boundary conditions while requiring that the 
diffeomorphism constraints hold on the boundaries of the spacetime. 

With these weak boundary conditions, the freedom in the action is vastly reduced as 
boundary terms like those of (|3.12j) must now be covariantly defined. That said, it is clear 
that the parameter space of possible actions remains large and, for example, includes all 
actions of the form 

I'h = Ih+ [ F s (^)+ / F B ( lij )+ [ F H (q AB ,u A ,l,e {e) ,v), (3.24) 

where the F are required to be dimensionally correct, but are otherwise freely defined vol- 
ume forms with the displayed functional dependencies. A lot has been written on possible 



boundary terms on B (see for example |ll|, |27j, |3Cj) but for now we simply note that the 



mathematical formalism allows these freedoms and that they can only be further reduced 
by imposing extra restrictions on the system. A good example of such a restriction is the 
requirement that the evaluated action not diverge in asymptotically flat spacetimes as B is 
taken to infinity. In this case extra boundary terms on B not only may, but in fact must be 
added to the action to remove such divergences. 

In this paper, the external boundary B is not the chief concern and so we will not consider 
these terms further. Instead we focus on the internal boundary H and so consider actions 
that are equivalent to In up to terms of the form 

AI = — J dvE[q, u, £_, 6(£),v] , (3.25) 
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FIG. 3: A three-dimensional volume along with its boundaries and their normals. 

where E[q, a), £ , 9^, v] is a free functional of the fixed data with dimension length. The 
requirement that the extra boundary term break up in this way is natural given the primary 
role played by the foliation of H . As we will now see, such terms are also very convenient 
as they allow for a clean Legendre transform of the action. 

E. Legendre transform to obtain a Hamiltonian 

In the previous sections it has been shown that the action In (|3.8|) is valid with either 
a strong or a weak fixing of the boundary conditions. In this section we introduce a notion 
of time to M which then allows us to perform a Legendre transform on the action and so 
derive a corresponding Hamiltonian functional. 

Flow of time and associated quantities 

Coordinate time must be externally imposed in general relativity. We do this in the following 
manner. First, foliate M with spacelike (partial Cauchy) surfaces {(£t, hij); ti <t< t 2 }; 
these will correspond to "instants of simultaneity" (see figure EJ). For simplicity, we assume 
that this foliation is compatible with all of the structures that we have already introduced. 
Thus E tl2 = Ex ,2 an d if {Bt',ti < t < t 2 } is the induced foliation of B, then B tl 2 = £>i >2 - 
On H we require that the induced foliation be compatible with S v and so it is natural to 
drop v and henceforth use the single parameter t to characterize both foliations. This extra 
structure on M and its boundaries allows us to extend the definitions of several quantities 
that were previously only defined on initial and final boundaries £> lj2 and Sx,2- Specifically, 
with the help of T, t we can define u a and n a everywhere on B and H. Further u a and 
7] := — u ■ n may be extended over B, while ( may be defined everywhere on H . 

In addition to a foliation parameter, we must also introduce a corresponding notion 
of evolution. This is provided by the evolution vector field T a which is assumed to be 
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compatible with the foliation in the sense that 



C T t = l. (3.26) 

Thus, the flow generated by T a evolves leaves of the foliation into each other in the required 
way. It is also assumed to be compatible with the boundaries so it maps them into them- 
selves; equivalently on B/H, T a G TB/H respectively. Then, on H it is natural to choose 
ilk) := T a - Note that nowhere in this construction have we assumed that T a be every- 
where timelike. Just as the natural (Killing) evolution vector field for the Kerr spacetime 
is not everywhere causal (even outside the event horizon), so in this case T a is allowed to 
be null or spacelike. Regardless of the nature of T a , it can be decomposed into components 
normal and tangent to the three-surfaces S t as 

T a = Nu a + V a , (3.27) 

where N is the lapse function and V a G TT, t is the shift vector field. Given that T a generates 
B and H we can perform similar decompositions on those surfaces. Respectively we have 

T a = Nu a + V a and (3.28) 
T a = = V a + V a , (3.29) 

where N = —T ■ u, V a = p/V A on B t and V a = p£V A on S t . Finally, the foliation 
parameter can be used to relate the various volume forms. Specifically, u a = —NV a t, so 
e = — iVdt A e. Similarly, on boundaries we have e = iVdt A e and e = dt A e, where 
dt is the appropriate pullback of dt. 



Legendre transform 

Using this foliation, we can break up the action into space and time components in the 
standard way jT^|. To this end, we introduce the momentum P ab conjugate to h a ^ as 

P ab = (Kf - K^h ab ) . (3.30) 

lo7rG 

Then, it is straightforward to show that 

NTZ = NH + V a H a + P ab C T h ab - 2C T Ku - K a h ab C T h ab + 2D a (h ab T c V c u b ) , (3.31) 



16vrG 

where 7i and 7i a are the usual Hamiltonian and diffeomorphism constraints, 

K = - 16nG (P ab P ab - P 2 /2) and (3.32) 

lb7rG 

H a = 2D b P a b (3.33) 

which vanish on shell. 

We can then substitute equation (J3.31|) into our expression for the action (J3.8|) and make 
use of the relationships between the various normals at B t and St given in (|3.1|) and (J3.2JI to 
express the action as: 



IH 



I -dt\H{t)+ [ P lj C T h ij + [ PC T e- [ PC T e\ . (3.34) 

J I J^t J St JBt J 
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Here the Hamiltonian functional H(t) is given by 



H(t) = I (S kNH + V*7U) + / { + W A u \ , (3.35) 



where 

p A a u b V a h b 



JA 8ttG 1 ; 



Furthermore, 

are identified as the conjugate momenta to hy, e, and e in the usual way (lH [l3L lisl Isil ] 



(E) - sinh 1 77 In C , 



T/ie Hamiltonian functional 

Although it is conventional practice to refer to H(t) as the Hamiltonian, it is clear that 
in many cases this functional is not really a Hamiltonian. The problem is that standard 
Hamiltonian mechanics properly applies only to closed systems. Here, however, energy flows 
are allowed through the boundaries B t and St and so in general H(t) ^ which in turn means 
that it cannot really be a Hamiltonian functional. Technically, the Legendre transform and 
Hamiltonian identification can only be correctly applied in those cases where the system is 
isolated — that is when H is an isolated horizon |Sj and B is similarly disconnected from its 
surroundings (a sufficient condition is that jC-Tjij — 0). However, in Hamilton- Jacobi 
arguments were used to demonstrate that the H(t) derived in this way can still be interpreted 
as the quasilocal energy associated with S 4 and its boundaries. These arguments are also 
supported by the extended phase space arguments of (lif . 

With this in mind, we can consider what H(t) tells us about the quasilocal energy of the 
system. First, since the bulk terms of f|3.35|) are pure constraints, the non- vanishing contri- 
butions to the energy come exclusively from the boundary terms. Thus, it is conventional 
to identify the value of each boundary contribution to H(t) as the energy of that boundary. 



In this case, the only such term appears on B t . This term was first identified in [ljj and has 
been much discussed since then. Note, however, that the value of the boundary energy is not 
uniquely specified. At B we are free to add any free functional of the boundary metric ^ a b 
to the action and hence obtain a different value for the energy of B. Indeed, in many cases 
(such as the limit of Bt going to infinity in asymptotically flat spacetime) it is necessary to 
add such reference terms so that the energy of B doesn't diverge. 

Next, we turn our attention to the horizon boundary. With our choice of action, we do not 
obtain a contribution to the Hamiltonian at H. However, as with the outer boundary B, we 
are free to add any functional of the fixed data to the action thereby obtaining a different 
value for the horizon energy. As we saw in Section ITnTjl this is a significant freedom. 
Specifically, if we impose the boundary conditions weakly, we can add any diffeomorphism 
invariant function of the fixed boundary data. The freedom (J3.25J) in the action gives rise 
to a corresponding freedom of 

AH = E[q AB ,u A ,Ji,d (t) ,v], (3.38) 

in H(t). Thus, as far as the mathematical formalism is concerned the contribution to the 
energy from St may be any functional of the fixed data. We will return to a more thorough 
discussion of the energy associated with a horizon boundary in Section El 
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First, however, we re-examine this problem from the canonical phase space point of view. 
This will both give us a new perspective on the energy calculation and also allow us to derive 
an expression for the angular momentum of trapping horizons. 



IV. CANONICAL PHASE SPACE FORMULATION 

In the previous section we have seen that trapping horizons can act as boundaries in a 
quasilocal action formulation of general relativity. Here we will see that they are also suitable 
boundaries in a canonical phase space description of gravity. Apart from the value that this 
demonstration has in its own right, it also gives us another point of view on some of the 
issues that arose in the last section. In particular, it will provide further elucidation of the 
relationship between the Hamiltonian and the energy functional H(t), and a sharpening of 
the discussion of the allowed freedom in that functional itself. Additionally, while the action 
calculation permitted us to investigate time evolution and hence obtain the quasilocal energy 
functional, this Hamiltonian investigation will allow us to consider more general evolutions 
of the boundary data and so also tell us about other "conserved" quantities. Specifically, we 
will be able to study angular momentum as conjugate to rotations of the boundary data. 

The particular formalism that we will use was described by the current authors in a 
recent paper and is designed to deal with situations, such as this one, where the energy 
of the quasilocal region is not expected to be conservered. Such systems cannot be dealt 
with using standard Hamiltonian methods which, by their very nature, are tailored to cases 
where energy is conserved. The methods that we will use are based on the generalized 
Hamiltonian formalisms of classical mechanics which are specifically designed to deal with 
situations like these where the energy isn't constant. These methods expand the phase space 
to include the time coordinate as a configuration variable (see for example Then, it is 

possible to construct a conserved Hamiltonian that generates time translations in this new, 
extended phase space. In this space, the energy of the system is no longer the value of this 
Hamiltonian, but instead is the value of the momentum that is canonically conjugate to 
time. In |l5( these methods from classical mechanics were extended to quasilocal general 
relativity. 8 Just as in classical mechanics, the time coordinate appears as a configuration 
variable and the energy as its conjugate momenta. However, for quasilocal general relativity, 
we also include coordinates parameterizing the two-boundaries in the configuration space. 
This allows us to determine the angular momenta of the boundaries, which are conjugate to 
rotations of the boundary data. 



A. Phase Space and Symplectic Structure 

The Hamiltonian approach to general relativity is based on a three-dimensional manifold 
E and its attendant fields which evolve in time with respect to some time parameter t, rather 
than the usual four-dimensional spacetime. In the case of interest here, E is diffeomorphic to 
[0, 1] x S 2 and has boundaries B and S (which are each diffeomorphic to S 2 ), as was shown 
in figure EJ Over this region, the basic defining fields of three-dimensional, time-dependent 
general relativity are the three- metric hij over E and its canonically conjugate momentum 



Other applications of extended phase space ideas to general relativity include the papers listed in |3, 
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P tJ := eP lJ . These variables contain all the information needed to resconstruct the full 
four-geometry at any time t and, once we specify a lapse function N and shift vector V\ 
also define how that geometry evolves in time (Appendix El . 

There are also degrees of freedom on B and S that determine how these boundaries evolve 
relative to S. To specify these evolutions it is sufficient to know the "orthogonal velocity" 
vector fields — that is v\- = (V • n)/N. These degrees of freedom are controlled by the 
area forms e and e and their conjugate momenta P and P which we defined in (IH.H7J1 . The 
encoding of this information is easily seen for B. There, 



N = and V i = {N^ri + V i , (4.1) 

and so v\- = r// ^/l + rf (recall that P is defined in terms of 77). On S the relation is a little 
more complicated, since there are two variables which combine to define the evolution; C 
relates horizon evolution to the null vectors and then £ ties those null vectors to £ and S. 
On S, 

N = a { C ~^) and ^ = a (c+^V + K (4.2) 

Then v\- = (2( 2 + C)/ (2( 2 — C) and the evolution in E is given by P (defined in terms of Q 
once we know C. 

In addition to these configuration and momentum variables, in the extended phase space 
we include the time foliation label t and spatial coordinates on both boundaries, 9 a and 9 a 
(for a = 1, 2), as configuration variables in order to allow for energy and angular momentum 
flows through those boundaries 0]. Along with these variables come their conjugate mo- 
menta Pt (a scalar field) and P Q and P a (one- form valued volume forms), which are closely 
related to energy and angular momentum. We can specify the evolution equations of the 
configuration variables; they are simply: 

£t = 1,4-^ = and 4-^ = 0. (4.3) 
dt dt dt v ' 

The values of the conjugate momenta and their evolution equations are initially left free, 
however we will see that the necessity of Hamiltonian evolution will greatly reduce that 
freedom — and so generate candidate expressions for the energy and angular momentum 
associated with the boundary B and horizon S. 

The allowed values of (hij, P u , e, P, e, P, t, P t , 9 a , P a , 9 a , P a ) span the phase space of 
Hamiltonian general relativity on the manifold M. The canonical symplectic structure 9 on 
this phase space then takes the form: 

n(S u 6 2 ) = J {(hP^fahij) - (5 1 P ij )(5 2 h ij )} + {(StPtWit) - (5iP t )(5 2 t)} (4.4) 
+ J {{8 2 P){8 x e) - (5 1 P)(5 2 e)} + J [{5 2 P){8{e) - (5 1 P)(5 2 e)} 
+ ^{(5 2 P Q )(5 1 ^) - (5 1 P a )(5 2 9 a )} + ^{(5 2 P Q )(5in - (5 1 P a )(5 2 9 a ) 



9 This is really the pre-symplectic structure since the existence of the Hamiltonian and diffcomorphism 
constraints means that we have not isolated the true degrees of freedom. The proper symplectic structure 
would restrict itself to the true degrees of freedom. That said, we will follow the standard practice and 
refer to this pre-symplectic structure as the symplectic structure. 
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Enforcing boundary conditions will impose restrictions on the allowed variations on S 
and B. There, we use the same weak boundary conditions that we considered in the action 
calculation in Section IIII Dl as well as imposing the appropriate diffeomorphism constraints. 
We now consider how the conditions translate into the canonical framework. 

First, on the outer boundary B, the three- metric 7^ is a combination of the lapse N, 
shift V A and two- metric Jab- We fix all of these quantities at the boundary, which in turn 
imposes restrictions on the full four-dimensional evolution at the boundary. From (|4.1j) we 
note that in particular this is sufficient to fix the bulk lapse N and shift V" 1 , up to the 
freedom encoded in 7] (or equivalently P). 

Next, we consider the inner, horizon, boundary S. The geometry of the horizon is nat- 
urally described in terms of the two-metric (Jab, shift vector V A , the quantities a and C 
and other quantities constructed from the null vector fields I and n. These can, of course, 
be related to the three-dimensional metric hij, the momentum P*- 7 , the lapse and shift, and 
the quantity £■ We will not give all the relations here, but refer the interested reader to 
Appendix [BJ In this case the boundary conditions are already tailored to the two-surface 
S. They require that we fix the two- metric qAB, the angular momentum one-form Co a and 
the value of aC (which is equivalent to the earlier fixing of £ ), as well as requiring that 0^ 
vanish. 

In analogy to the way that we restricted the variation of these quantities in the action 
formulation, here the variation of "fixed" quantities is restricted so that 5 respectively acts 
as: 

5xs = i ~ 5t) It +Crx and ^ = (^ + £«r, ( 4 - 5 ) 

on S and B. Here, St is a function of t alone (which means that the variation preserves the 
foliation) and 5X/5X G TB/S. Finally, note the difference between this decomposition of 
the variation and that given in (|4.5|) . In the earlier version, we decomposed with respect to 

. a ~ 

V rather than T. Thus, the two are related by: x Q = St and 5X = x A — (5t)V A . 



B. Hamiltonian evolution 

In this section, we will use Hamiltonian methods to obtain expressions for the energy and 
angular momentum of a horizon. In order to do this, we will be interested in evolutions of 
the form 

5 A = X j t + C x , (4.6) 

where A Q = A„(£) is a scalar function and A' G TS is a vector field defined over E. In 
addition, we require that A* be tangent on the boundaries. That is 

X i =p2~X A onS and A* = pj\ A on B , (4.7) 

for some X A G TS and X A G TB. For an evolution along A, we will show that there exists 
a Hamiltonian functional which generates the evolution. Specifically, this means that the 
symplectic structure (|4.4|) with Si = 5\ can be manipulated into the form 

£1(5 a, 5) = SKa + constraints . (4.8) 
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This class of evolutions includes both time evolutions and rotations of the system, so we can 
expect the Hamiltonian functional K A to provide information about both the energy and 
the angular momentum associated with the horizon S, subject to an appropriate choice of 
A. 

In order to rewrite the symplectic structure in the form (J4.8)) we must evaluate the bulk 
term by making use of the evolution equations for the three-metric and extrinsic curvature 
derived from the Einstein equations. This is a long calculation and we will simply state 
the result; we refer the interested reader to Appendix El for more details. Before giving the 
result, it is useful to introduce a shorthand 

H A = f {X (NH + V^Hi) + X'Tii} , (4.9) 

for the bulk part of the Hamiltonian, where 7i = eH and Tii = eHi. As usual this is 
comprised of constraints and will therefore vanish on shell. Then, we can express the bulk 
contribution to the symplectic structure as 



L 



{(5P^')(5 A %) - ((J A P«)((5/ly)} (4-10) 

= S (* A - L {^G mK + e{X ° VA + XA)U ) + 8^ \ 
+ J { (S A P) (Se) - (5P) (6 A e) } + J { (<J A e) (6P) - (SI) (6 A P) } 

+ ^ L{^ +euvA } +m L c 4^G ekn+euvA . 

+ J {(5X a )ej a + (S9 a )S A (ej a )} - J {(5X a )eCu a + (69 a )6 A (eu a )} . 

In order to obtain the above, we have made repeated use of the boundary conditions. Specif- 
ically, we have used 9m = and the fact that all other "fixed" quantities vary as described 
in (|4.5|) . Additionally, we have made use of the diffeomorphism constraints on both the 
horizon and outer boundary. Then, ()4.10|) is close to the form that we require. The first 
term is an exact variation, and will be a contribution to the Hamiltonian K A , while the 
second line will vanish provided the bulk constraints are satisfied. The third line consists 
of terms involving the boundary area forms and their conjugate momenta. The final two 
lines involve additional terms on the boundaries. These terms will eventually lead to our 
expressions for the energy and angular momentum of the horizon and outer boundary. 
To shorten and clarify our expressions we introduce further notation. First, we define 

K t :=P t -[e l-^Nkn + V A j A \ , (4.11) 



JB 

L A ■= Pa - ej A and L A ■= Pa + ^— ■ 

07TG 
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We then substitute ()4.1(J|) into the symplectic structure (|4.4|) . Then, 



n(6 A , 6) = 5\H A + \ D K t + J \ A L A + J ~\ A ~L^j (4.12) 

- / d 3 x {s(\ N)n + 5(\ v i + x i )n i ] 

-(SX )K t - [ (SX a )L a - / (SX a )L a 
Jb Js 

-(5t)(5 A K t )- [ (59 a )(5 A L a )- [ (59 a )(5 A L a ) . 
Jb Js 

The first term line in ()4.12j) is an exact variation. Therefore, if we can argue that the 
remaining terms vanish, this term will be the Hamiltonina K A which generates evolutions 
along A. This will be the case provided the five quantities T-C,7ii, K t , L a and L a are all 
zero. In addition, the evolution along A of K t , L a and L a must also vanish. 

The vanishing of 7i and Tti is nothing more than the usual Hamiltonian and diffeomor- 
phism constraints. The other conditions restrict the values of the momenta Pt and P a . Until 
this point, these quantities have been entirely free, but now we must require 

(413) 

P Q = ej a and V a = -——u) a . (4.14) 

07TG 

Next, the fact that the A variation of K+ and L n vanishes fixes the evolution of the momenta. 



Specifically, it guarantees that the functional forms given in 1)4.13)1 and 1)4.14)1 are preserved 
in time. Thus, our evolution will be generated by a Hamiltonian provided that the standard 
bulk constraints of general relativity are satisfied and in addition the conjugate momenta 
introduced on the boundaries take the form given above in ()4.13j) and ()4.14|) . 

Thus, if these equations of motion are imposed in the extended phase space, 5 A evolutions 
are generated by the Hamiltonian 



A a = H A + X K t + / e\ A L A + / e\ A L A . (4.15) 
Jb Js 

As always for a Hamiltonian, this is constant on-shell. In fact it vanishes since it is composed 
entirely of constraints. However, in the extended phase space, this does not mean that the 
energy and angular momenta are identically zero. Instead, the energy is given by the value 
of the momentum canonically conjugate to the time variable. Thus, the formalism tells us 
that the energy associated with the evolution is 

The value of energy we have obtained from the phase space calculation agrees with the one 
arising from the Legendre transform in the previous section. Again, we see that there is no 
horizon contribution. Similarly, the angular momenta of the boundaries are determined by 
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the momenta P a . Therefore, the angular momentum J[0] associated with a vector field <fi A 
at the two boundaries is 

J B [4>] = [ e4> A j A and J s [cf>} = --L / e^W (4.17) 
Jb 8?rG Js 

We thus obtain angular momentum values for both boundaries independently. The value at 
the outer boundary agrees with what is generally obtained. The expression for the angular 
momentum of the horizon provides an extension of the previous definition (J2.6j) on an isolated 
horizon. This phase space calculation tells us that we can associate an angular momentum 
with any vector field <fi at the horizon, not just axisymmetric Killing vectors. As expected, 
the value of the angular momentum is determined by the field uja- 

Freedom in the Hamiltonian 

Just as there was a freedom in the definition of the action Ih and the quasilocal energy H(t), 
there is a corresponding freedom in the definition of our Hamiltonian K\. A close examina- 
tion of (|4.12|) reveals that K\ is not the only Hamiltonian that generates the evolution <5a 
through the physical phase space. Consider for example 

K' A = K A + \ F(t), (4.18) 

where F(t) is a free function. Then, it is fairly easy to see that this functional also generates 
the usual physical evolutions subject to a modification of the constraint equation ()4.13|) so 
that it becomes 

p > = IA^ +9A,A } +F(t) - (419) 

In the extended phase space, this corresponds to a displacement to an alternate constraint 
surface that still corresponds to the usual physical evolution of the metric and extrinsic 
curvature. In fact one can argue that this is the most general term that can be added to 
the Hamiltonian ^5[. With this alternate Hamiltonian, the quasilocal energy becomes 

E[d/dt] = F(t) + J^{^§ + V A 1a} ■ (4.20) 

The quasilocal angular momentum remains unchanged. 

Mathematically, F(t) can take any form but from a physical point of view the only 
functionals that make sense are those constructed from from the time dependent boundary 
data. That is 

F(t) = Ffr AB , N, V- q AB ,Co A , aC, 9 {£) ] , (4.21) 

and so the quasilocal energy derived by this method is compatible with that ()3.38|) found 
by action arguments. 
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V. ENERGY AND ANGULAR MOMENTUM OF THE HORIZON 



In Sections IIHI and IIVI we have obtained expressions for the energy and angular momen- 
tum of a spacetime containing a horizon as a boundary. In this section we will consider 
these expressions in more detail. In particular we will show that the expression for angular 
momentum is equivalent to the corresponding expressions on isolated and dynamical hori- 
zons. Turning to the energy we will see that although a wide range of energy expressions 
are allowed by our mathematical formalism, this freedom may be greatly reduced through 
the imposition of a few natural conditions. Once this is done, the remaining freedom is 
essentially that found in the dynamical horizon formalism. 



A. Angular Momentum 

The angular momentum of the horizon, derived using the Hamiltonian framework, is 
given by 

This expression is unique, as the freedom we have to add terms to the Hamiltonian K\ will 
not affect the value of the angular momentum obtained. Furthermore, this value of angular 
momentum is applicable for both isolated and dynamical horizons. It is immediate from 
(12. 6 j) that our expression is identical to the standard one on an isolated horizon. 
The dynamical horizon definition for angular momentum (|2.7j) is 

Jdh[4>] = 77^77 / e^f'Vah. (5.2) 
87tG J s 

This is clearly very similar to the expression that we derived by Hamiltonian methods, but it 
isn't identical — (J5.1|) made use of the null normals l a and n a whereas ()5.2|) uses the spacelike 
normals f a and f a . We will now show that for reasonable choices of (fi A , this difference is 
irrelevant and either choice of normals will give the same value for the angular momentum 
(up to an overall sign) associated with <p A on the horizon. 

To see this, let us consider the conditions that <p A should meet so that it may be viewed 
as a generator of rotations. We start with the example of a two-dimensional slice of the Kerr 
horizon. There, the obvious choice for (fi A is the Killing vector field, (<9/<90) in the usual 
coordinate systems. Note that in addition to satisfying C (d/d^QAB = 0, this vector field can 
also be integrated to foliate the two-slice into a family of closed curves, all of which have a 
common period. 

With this example in mind we will assert that <p A is a suitable "rotation" vector field if : 

1. It integrates to a congruence of closed curves that foliate S (rotations should be peri- 
odic!). 

2. C t^e = (Ia4> a = 0. Given the first condition, this one may always be met through 
a rescaling of 4> A . Specifically, (fr A is restricted to take the form 

<p A = pi AB d B e , (5.3) 

where 9 is a labelling of the congruence and — (3(6) is a free function. This is a 
necessary condition for (f> A to be a Killing vector field. It may be thought of a requiring 
a kind of uniformity of <p A along each curve. 
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3. Finally, (3 should be chosen so that the curves all integrate to have a common period 
(rotations should have a common period for all orbits!). 

A cf) A of this type satisfies many, but obviously not all, of the properties of a Killing vector 
that generates rotational symmetries. It is not hard to show that these conditions are more 
than sufficient to force the two angular momentum expressions to be equivalent. Begin by 
noting that on a (spacelike) dynamical horizon, 

^r(f a + r a ) and n a = ^=(f a - r a ) . (5.4) 



2 v ' ^2C 

Whence, it easily follows that 



s 



J~e^ a hV a f b + l -d A { ( f ) A \nC) - hnCd A ^ 



The second term on the second line vanishes because S is closed and the third term vanishes 
since <p A is required to be divergence-free. 

Thus, with <f> A of the class that we have discussed, the angular momentum expression that 
we derived by Hamiltonian arguments ()5.1|) agrees exactly with the expression that appears 
in the Ashtekar-Krishnan flux law ()5.2|) . The two calculations then support each other. 
Our derivation gives a Hamiltonian justification for viewing this expression as representing 
angular momentum, while the flux law then shows how to calculate changes in that angular 
momentum due to matter and energy crossing the horizon. 

There are a couple of points to keep in mind with respect to these expressions. First, 
the derivation of the dynamical horizon flux law does not actually make use of the fact 
that if is a dynamical horizon, except for assuming it is spacelike (so that r a and f a are 
defined). Thus, it is true for all spacelike surfaces, not just horizons. Second, although [j| 
is the first paper to derive a flux law of this type for spacelike surfaces, very similar angular 
momentum definitions and flux laws have appeared for timelike and null surfaces. In fact, 
while expressions for quasilocal energy vary widely, angular momentum is almost always 
defined with respect to the connection on the normal bundle to a two-surface. Examples of 
these expressions include: the Komar formula for stationary spacetimes (see, for example, 
j26|), the definition of angular momentum for an isolated horizon and the Brown- York 



definition for quasilocal angular momentum The only difference between the various 



definitions is which normal vectors to S are used in the angular momentum expression. 



B. Energy 

Both the action and Hamiltonian calculations allow for the horizon energy to be any free 
functional of the fixed data. That is, 

Energy = E[q AB , Cu A , 6> w , £_, v] , (5.6) 

where the exact form of E is unspecified. However, it is also clear that many of these 
functionals will not be physically acceptable. To narrow down the choices, we impose the 
following conditions on the energy functional: 
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1. Dimensional correctness: When c = G = 1, E should have units of length. 



2. Constancy in isolation: E should be constant over isolated sections of horizon. 

3. Scaling with time: E should rescale linearly with reparameterizations of time. 

All three of these are very natural. In particular any functional that doesn't satisfy the first 
two conditions could hardly be considered an energy. The third is a little less precise, but 
is also very fundamental. Energy is conjugate to time and so its value should change with 
rescalings of time. However, what has been left imprecise in this statement is exactly what 
we mean by time. We will now consider this point. 

On the outer boundary B, the evolution vector field T a is naturally identified with the 
flow of time. There, T a appears in the energy functional 

Nku , -TrA- \ f ~ ma I ku^a 



E ^JA-^ +vu ) = L~ eT '{-i^ +s °)- (57) 

Clearly this is dimensionally correct and also rescales with T a . Further, since T a is fixed data, 
any extra reference terms that are added may also be given the correct scaling behaviour. 
Given this experience, on H the obvious first candidate to consider as a time vector is (d/dv) 
or possibly V. However, with the horizon boundary conditions, this cannot be the correct 
choice as these vectors neither appear in our base expression (which is on H) nor in the 
fixed data. The only way that we could get a rescaling with respect to these would be 
through the foliation parameter v. However, an expression of this type would clearly fall 
afoul of condition 2 and not be constant on isolated horizons. 

Thus, the energy associated with the surface S cannot rescale with either (d/dv) or V. 
However, among the fixed data we do have £ = aCdv. We are then naturally lead to 
consider an energy functional that has a linear dependence on £_. Furthermore, the null 
vector field i a is causal which gives it an additional attraction over other, possibly acausal, 
alternatives. In considering possible forms of the energy functional, we note that over isolated 
regions where i a is normal to H, £ = 0. Therefore, we propose an energy functional of the 
following form: 

1 



dE = — / £ Ae , (5.8) 
c a a Js *~ 

where a is the area of S and c Q is a constant. That is, we argue that the energy flux should 
be proportional to the "average" value of £ . 10 Then, given cross-sections Si and S2 which 
bound a region AH their energies are related by: 



1 r . ji 1 



E 2 -E l =\dv{ — a? = - / aCdv, (5.9) 

. c o a J S J c o 



('2 



where aC = i f s 'eaC is the average value of aC over any surface. Thus, for any choice of £, 
an energy defined in this way will meet our minimum requirements: it will be dimensionally 
correct, constant when H is isolated, and (its flux) will rescale with respect to a causal 



Note that our freedoms would equally well allow us to multiply the right-hand side of this equation by 
any dimensionless functional f[qAB, v] of the remaining fixed data. However, since t may already be 
freely chosen, this does not add to the generality of the expression. 
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"time" vector field. However, the energy in ()5.9|) does not fit naturally into the allowed 
freedoms of the action formulation — we would have to define E as an integral of the fixed 
data from some reference surface S a with energy E Q to the surface S — and is even less 
natural in the Hamiltonian framework. Therefore, we will further restrict the freedom in £ , 
and hence the form of the energy. 

For a specific horizon geometry, we would like to fix a unique normalization for I and 
consequently a preferred value for the energy. This is familiar from, for example, the Kerr 
isolated horizon. There we obtain a family of permissible evolution vector fields £, but 
restrict to a preferred one by requiring it to have a specific value of surface gravity. To make 
a similar fix here, we need to tie £_ to the geometry of the horizon. If the horizon does not 

contain any "partially-isolated" two-surfaces 11 then a nice way to do this is to require that 
£_ satisfy 

£_ = c o df[q AB ,0JA] aC = c D ^jQ, (5.10) 

where / = /[^abj&a] is some functional of the fixed geometric data on each two-surface 
S. Recently it has been shown that a wide class of 6^ = three-surfaces actually meet 
this assumption j35|. Roughly speaking, if there are trapped surfaces arbitrarily close to 
(though spacelike-separated from) each two-surface then on those {6m = 0) surfaces £ 
vanishes either everywhere or nowhere. Thus, for most H which can be identified as black 
hole boundaries, one should be able to fix £_ in this way. 

With such a fix, £ is "constant" over each leaf of the foliation. Then, consulting (J5.8|) it 
is not hard to see that this will uniquely determine the flux of energy through the horizon 
between cross-sections - - and so the actual energy expression up to the usual additive 
constant. However, in general relativity the only dimensionful constants are the speed of 
light c and Newton's constant G. From these, it is impossible to construct a quantity with 
units of energy whence it is natural to set the additive constant to zero. Then, we obtain 

E = f[q ABl u A }. (5.11) 

Thus, with this choice the energy is fixed directly in terms of the geometry and £_ is fixed 
up to an overall constant rescaling. 

The energy expression derived above has another interesting property — it must always 
be non-decreasing. This follows immediately from the fact that the pull-back of the forward- 
in-time pointing I will always be positive on a dynamical horizon H (or zero in isolation). 
Then, with Kerr black holes and the Penrose process in mind, this suggests that we associate 
E with the irreducible mass of the horizon (see for example rather than the "total" 
mass which can either increase or decrease in time. That is we define 

/ = (5-12) 

where R is the areal radius defined by a =: AirR 2 . 

If one wished to instead consider the "total" energy it would be necessary to consider a 
more complicated evolution vector field. Specifically in the case of axisymmetry, it would 
be natural to consider I + where is the rotational Killing vector field and Q is a 



That is, two-surfaces which contain both regions where £ = and others where £ ^ 0. 
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function of the fixed data. Then the "irreducible" energy would remain unchanged since 
t + Q(f) A e = ^Ae. However, <p A would enter the fixed data as a symmetry of the 

two-metric, and so one could also include it in energy expressions. In particular, with this 
symmetry there would be an unambiguous angular momentum J for each cross-section of 
the horizon and this could be used to construct energy expressions like the usual Smarr 
formula E(a, J) for the total mass of a Kerr black hole. 

Relation to dynamical horizons 

We will now show that the energy obtained by Hamiltonian methods is related to the 
dynamical horizon energy and flux laws. Furthermore, by appealing to dynamical horizons, 
we can fix the unknown constant c Q and hence further restrict the energy. We begin by 
reviewing those aspects of the dynamical horizon calculation which are relevant here. 

Dynamical horizons are required to be spacelike, therefore, we can introduce the unit 
timelike and spacelike normals f a and f a . From these, we obtain a null vector 

t = f a + r a . (5.13) 

The dynamical horizon flux law is obtained by integrating a combination of the Hamiltonian 
and diffeomorphism constraints over the horizon. The specific combination is chosen to 
correspond to the £ a f b component of the Einstein equations. Then, by making use of the 
fact that at the horizon 9(g) = 0, one obtains an energy flux law 

E 2 ~ Ex = J iH eT ab C E f b + ^ J iH) eN E {af e fa mB + 2( A ( A } . (5.14) 

In the above, {H) e is the horizon volume form, (a = P^^V^a, N E is the lapse on the horizon 
and £ E is the evolution vector field associated with the energy functional E. Finally, the 
energy and evolution vector are tied together via the lapse by 

N E = 2\\dE\\ and £ a E = N E t. (5.15) 

The expression (|5.14|) provides a nice relationship between the change in energy of a horizon 
and the flux of gravitational and matter energy through the horizon. The change in energy 
is shown to be equal to three terms, all of which are non-negative (provided the weak energy 
condition is satisfied). The first is identified as the flux of matter through the horizon and 
the second is the flux of gravitational shear. There is no clear interpretation of the third 
contribution. 

Thus, the dynamical horizon energy always increases in the same way as the "irreducible" 
energy that we have derived. Further, we note that ()5.15|) can be rewritten without reference 
to unit vectors (with the help of (|5.13|) ) as 

£ =2dE. (5.16) 

This suggests that we fix the unknown constant arising in (|5.1U|) as c Q = 1/2. Then, the 
Hamiltonian evolution vector field l a is equal to the evolution vector £ a arising in the flux 
law f)5.14j) . Therefore, the Hamiltonian evolution vector field and energy also satisfy the 
dynamical horizon flux law. 
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In summary, a simple argument has lead us to obtain an energy functional for the horizon. 
This energy is constant whenever the horizon is isolated. Additionally, by arguing that the 
energy should rescale appropriately with reparametrizations of time, we have identified the 
null vector £ as the evolution vector field of the horizon. This is a surprising result. One 
would usually expect the evolution vector field to be tangent to the boundary, but in this 
case it will not be, unless the horizon is isolated. Next, by relating the value of £ to the 
other geometrical properties of the horizon, we obtain a energy associated to £ which is 
unique up to an overall constant rescaling. It is non-decreasing and so naturally associated 
with the irreducible (rather than total) mass of the horizon. Finally, by appealing to the 
dynamical horizon results (which also describe non- decreasing energies), we can restrict the 
little remaining freedom in the energy in such a way that £ and Ef satisfy the dynamical 
horizon flux laws. 

VI. DISCUSSION 

There are many different local definitions for black hole horizons: apparent, trapping, 
dynamical, isolated, and slowly evolving. Although not identical, these all share some com- 
mon features. Specifically, they are all three-surfaces which can be foliated by two-surfaces 
of zero outward expansion. More precisely, for each two-surface the outward directed null 
normal £ to the surface has zero expansion, 9u\ = 0. In this paper, we presented both 
action and canonical phase space formulations for a spacetime with such a boundary. These 
formulations provide us with quasi-local expressions for the horizon energy and angular 
momentum. 

A trapping horizon has one natural boundary condition, namely that the outward ex- 
pansion 9(t) vanishes. This in itself is a mixed boundary condition in the sense that it 
depends upon both the intrinsic and extrinsic geometry of the surface. We impose addi- 
tional boundary conditions which fix the two-metric <Jab and the one-form uja on every slice 
of the horizon. The one-form uja is the connection on the normal bundle to the horizon. 
Physically, these fix the intrinsic geometry of the two-surfaces as well as their associated 
angular momenta. Unlike the 8^ = condition which holds for all rescalings of the null 
vector £, the one-form ua depends upon the particular normalization of £. Thus, we must 
also fix a preferred normalization for £, and this is done by fixing the pull-back of £ to the 
horizon. 

We have shown that these are suitable boundary conditions for both the action and 
Hamiltonian formulations of general relativity. While this is an interesting result in its own 
right, in this paper we are mostly interested in the expressions that we can then derive for 
the energy and angular momentum associated with the horizon two-surfaces. The result 
for the angular momentum is unambiguous: there is a definite value of angular momentum 
associated with each cross section of the horizon. Furthermore, this value is identical to the 
angular momentum calculated for axisymmetric isolated horizons. In addition, if we restrict 
attention to divergence free vector fields, the value of the angular momentum agrees with 
that arising in the dynamical horizon flux law. 

Obtaining an energy for the horizon is not so straightforward. As in all quasilocal Hamil- 
tonian calculations, there is a freedom to add a functional of the fixed boundary data to 
the Hamiltonian and hence change the value of the energy. However, we can almost entirely 
remove this freedom by requiring that the energy satisfy a few simple physical properties. 
These then lead to several interesting results. The first is that the appropriate "evolution 
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vector" for the horizon is the null vector i and hence, in contrast to a typical Hamiltonian 
formulation, the evolution vector need not map cross sections of the horizon into each other. 
Next, by requiring that the null vector £, or more accurately its pull-back to the horizon, 
is compatible with the horizon geometry we obtain a range of possible expressions for the 
energy of the horizon: For a particular choice of null vector field i the horizon energy is fixed 
up to an overall constant rescaling. All of these energies are constant on isolated horizons 
and increasing on dynamical horizons, and can therefore be interpreted as describing the 
irreducible mass of the horizon. Finally, appealing to the dynamical horizon formulation we 
can then fix the value of the rescaling constant relating the energy to the evolution vector 
£. Once this is done, we find that our energy Ei and evolution vector field i satisfy the 
dynamical horizon flux law. 

Note however that this is certainly not meant to imply that all problems concerning 
energy and dynamical horizons have been resolved. While it is remarkable that the allowed 
energy expressions obtained by these two very different methods agree, it is also important 
to keep in mind that neither method singles out a unique expression. While this is somewhat 
disappointing it is also not too surprising given the well-known ambiguities of gravitational 
energy. 

Other problems also remain for future investigations. Probably the most obvious arises on 
recalling the discussion of footnote El There it was pointed out that these 9^ = horizons 
are not uniquely defined — in most spacetimes a given horizon will be part of an entire 
family of 9^ = three-surfaces which may be smoothly deformed into each other. Thus, 
in using these techniques to discuss energy flows, one needs to be careful as quantities such 
as cross-sectional area will certainly depend on the member of the family chosen. This is 
not necessarily a problem. In physics, it is common for observations to vary with reference 
frames (think special relativity). The important thing is that invariant physical predictions 
should be compatible. Thus, for example, while the total energy absorbed by a black hole 
during a finite incidence of gravitational waves should be invariant, the details of how that 
happens may depend on the specific "horizon" . 

Lastly, the formalism and boundary conditions presented here may find other applica- 
tions. Perhaps the most interesting is the possibility that the boundary conditions that 
we have discussed could find a use in quantum gravity discussions of dynamic black holes 
- just as the isolated horizon boundary conditions and Hamiltonian formalism found an 
application in the loop quantum gravity computations of black hole entropy. 
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APPENDIX A: DETAILS OF THE ACTION VARIATION 



In this section we show that 



SI H = — *— [eG ab 5g ab l — I ( % (K% - K a h ij ) Shu (Al) 



+ 



b S? « - K ^ ^ + /<*> 8inlrl <") 



+T7T^ / d«Ae{^} + /(5e) (1 - InC) 

as asserted at equation (jXTUJ). In the interests of brevity, the (well-known) outer boundary 
terms on B will be omitted from this derivation. Details of the calculations for those terms 
can be found in a standard reference such as 



ism 

We begin by writing the bulk term in the action as: 

5 {ell) = e {G ab 5g ab + V a [g bc (V b 5g ac - V a 5g bc )] } (A2) 

On integrating over M and using Stokes theorem to move total derivatives out to the bound- 
aries this becomes: 

r ell) = [ eG ab 5g ab + [ iJ:) eu a g bc (V b Sg ac -V a 5g bc ) (A3) 

AI J JM JY, 2 -Y,i 



dv A er a g bc {V b 5g ac - V a Sg bc ) 



H 



In order to obtain the term on the horizon, we have used the expression (|3.6|) for the volume 
form. With just a little algebra, the term on S 12 can be written as: 

[ {S eu a g bc (W b 5g ac -W a Sg bc ) (A4) 

= - / { { %{K i i-K^)8h ij + 28C ) eK il )\- [ eh a 5u a ■ 

JS 2 -^i ) JS2-S! 

The last term is not, of course, zero since variations of the "time-space" components of g ab 
can give 5u a a component in the n a direction. This is in contrast to Su a whose "direction" , 
though not normalization, is fixed by the requirement that it be normal to the S 12 surfaces. 

The signature on H is indeterminate, so in analyzing the terms in (|A3|) on that boundary 
we must avoid the use of the inverse metric as well as proper volume forms or unit- 
normalized normal vector /one- form fields. With these caveats in mind we note that: 

r a g bc (V b Sg ac - V a 5g bc ) = (r a q bc + n a tV c ) (V b 5g ac - V a 5g bc ) (A5) 

In turn, the two expressions on the right-hand side of this equations can be rewritten as: 

r a q bc (V b 5g ac - V a 5g bc ) (A6) 
= -25k T - 5q ab K b + q ab V a 5r b + q b a V a 5r b - {5q ab )q ac {tn d + n b £ d ) V c r d , 
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where k T = g ab V a Tfe, and 

n a £ b V c {V b 5g ac - V a 5g bc ) (A7) 
= -25k v + 2(5V>i - V\5n a V b t + n a V b 5t + 5n a V b £ a + n a V b 5£ a ) . 

There is still significant work needed to transform these expressions into what is present in 
fjAljl . We begin by considering the right hand side of equation (IA6J) . First, we note that the 
5k T term can be rewritten using: 

ek T = eq ab V a r b = -eq ab V a V b + 2ea9 {£) = -£ v e + 2ea9 {i) (A8) 

Then the L ye can be integrated out to the future and past boundaries of the horizon while 
the second term (and its variation) will vanish since we require that 9^ = on the horizon. 
The second term in (jA6|) appears as a contribution to s AB 5qAB in the final expression 12 . 
Next, we turn to the third and fourth terms in equation (jA6|) . those containing a 5t. By 
expanding r in terms of £ and n using ()3.5|) , it follows that they can be rewritten as 

q ab V a Sn + ~q a b V a 5r b (A9) 
= 2(5a)9 (i) + 25(aC)6 [n) + g afe V Q (g bc 5r c ) - (q ab C v q ab )(n a 5£ a + n a 8t) • 

The first two terms arising above appear in the final formula (jAlJ) and the third is exact and 
will therefore vanish upon integration over a two surface S v . The final term will combine 
with others to form an exact derivative. To see this, we add the last term on the right hand 
side of (jA6|) to the last term of ()A7|) . By making use of the fact that n a V b T a = n a V b V a and 
£ a V b r a = -£ a V b V a it follows that 

-(5q ab )q ac (£ b n d + n b £ d )V c r d - V b (5n a V b £ a + n a V b 5£ a + 5n a V b £ a + n a V b 5£ a ) 

= -C v (n a 5£ a + n a 5t) . (A10) 

When added to the £ y term from (|A9J) we obtain an exact derivative which can be integrated 
out to the future and past horizon surfaces Si and S2. 

Finally, we turn our attention to the expressions in (jA7|) . The first appears in the s AB 5q~AB 
of the final expression (through q AB 5q~AB = (j ab ^Qab) while the last was dealt with in (jA10|) . 
Therefore, we are left only with the 2{5V l )uJi. This can be rewritten using 5(V l dv) A e = 
— Pa (8V A )dv A e to show that 

2 dv A e(5V l )uJi = -2 5 (dv A eV A u A \ - 2(5 dv) Aeny + 2 V A 5(dv A eu A ) . (All) 

The results from fIXsl- IXTTT) can be used to rewrite the horizon terms from (IA3I) as: 

Jdv Mr a g bc (V b (5g ac ) -V a (5g bc )) = -2 5 (jf {du A e (k v + V a Cj a ^ + 2 nd (i) A e}j (I 

+ [ dvAe [(q AB k T - k AB + q AB K V ) 5q AB ] 

JH 

+ J ^2(5 V) Ae6 [n) -2(5n) Ae9 {i) + V A 5(dv AeCo 



+ f {2 5I-e[n b 5£ b + n b 5£ b ]} . 



12 To be precise, note that k^ b 8q ab = -k Tab 5q ab = -k Tab 8{p^p^q AB ) = -k rAB 5q AB = k AB Sq A B- 
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In the above, = Pi a £ a = aC[dv]i and (n)i = Pi a n a = a[dv]i. By using the expression 
(I3.11J1 for s AB we see that contribution on H is of the form arising in (jAljl . leaving us only 
the corner terms arising at S± and S2 to consider. We have a contribution rit,5£ b + n b 5£b 
arising from H and another contribution, h a 5u a from S 12 . However, there are relationships 
between the unit vectors u, h and £, n given in (|3.2|) . Making use of these, and the fact that 
5u a oc u a , it follows that 

h a 5u a + njt + n a <% = -25(ln C) • (A13) 

Then, combining (|A12|) . ()A13|) and recalling terms in the action (|3.8|) on £ lj2 and H, it is 
straightforward to arrive at the desired result, (jAlj) . 



APPENDIX B: MANIPULATING THE BULK SYMPLECTIC STRUCTURE 

In this appendix, we will rewrite the bulk symplectic structure in a form which is directly 
useful in the body of the paper. We are interested in the case where Si = 5 A generates 
an evolution along the vector field A a . In three-dimensional language, the evolution of the 
three-metric, hij and its conjugate momentum P* J which encodes the extrinsic curvature 
degrees of freedom are given in by: 

hhi^Kjhij + Cxhij and 6 A P^ = A ^P iJ + £ X P^ . (Bl) 

If the bulk fields are required to satisfy the Einstein equations, then we can easily derive 
expressions for the time evolution of and P* J as: 

jhij = teKGN{Pij - X -Ph l3 ) + £ vKj . (B2) 

_P*i = — ( N(R ij - -Rh lj ) - (rrrPN - h ij D k D k N) ) + £ v P ij (B3) 
at IQttG \ 2 / 

+8nGN ((P kl P k i - ^PP)h ij - 4(P fc(i P fc i} - ^PP ij )^j , 

where R is the three-dimensional Ricci scalar and P = hijP^ and N and V a are the lapse 
and shift respectively. Then, a fairly standard calculation 0, HH allows us to rewrite the 
bulk symplectic structure as: 

/ {(5P^)(5 A h tJ ) - (5 A P^)(5h tJ )} (B4) 
= J {\ N5H + (\ l +X V i )6H i } 

l 5q AB 

+ [ 2S(eh i P ijP J)(\ V A +~\ A )+ [ \2X (h i V i )5(eP^h i n j ) - ^5{ek h ) \ , 
Jas JdT, I 8ttG J 

where in the interests on conciseness we use tildes to represent two-dimensional quantities on 
all boundaries. Thus, Pab = PaPb^j * s the pull-back of P^, V A and X A are the components 



\ (£> I/^P n CA n DB I (h a AB h AB \ I ^o^ %D i N ~AB 

- I e\-\ [n k V )P C Dq q + ~[QtiG ~ 16ttG 
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of V 1 and A* in T(<9E), and qAB is the two-dimensional metric on <9£. Additionally, since n 
is a unit vector in E but orthogonal to 9E, k AB is the extrinsic curvature of the boundary 
<9£ in E. 

Now, we would like to rewrite the boundary contributions in terms of quantities which 
are tailored to the boundaries B and H. In the four-dimensional context, we are more 
interested in the normals to B and H in M rather than n, the normal to <9E in £. The 
former are independent of our spacetime slicing while the latter depends upon the foliation of 
the spacetime. To facilitate this procedure, we re-express the familiar boundary quantities 
from Section IIHI in terms of the three-dimensional quantites appearing above. We begin 
by writing out the various horizon quantities in terms of the three-dimensional metric and 
extrinsic curvature as: 



0(n) 



.AB 



~AB 



i.l 



n n 



P 2 /2] + WDiN 



(WirG^jPijti + d A (lnC) 
^(lnC) - CylnC + 16nG(V k n k ) [P< 
C [-(l&irG)P ij h i n i + Dih*] 

(lGirGn^q^fvpcD + Nk AB 

-(16irGV%)q AC q BD P D + N {k h q AB - k AB ) + {fjDiN)^ 3 



+q AB 



|(lnC)-£v(lnC) 



Similarly, we can write the outer boundary terms as: 

k AB = y/l^rfk AB + (WnGri) \i ac i bd P C d ~ P^ AB /2 



(B5) 
(B6) 

(B7) 
(B8) 

(B9) 
(BIO) 



n c A c N 



3a 
Ns AB 



n 



DiN + l6irG(V k n k ) [n'tiPPy - P/2] + — (sinh" 1 //) - C ^smli" 1 r/)(Bll) 

(B12) 



— IQnGp^Pijhi — G?A(sinh x ?7) 
—(167rGfiiV i )j AC j BD PcD + iV(^7 AB - O + {fjDiN)^ AB 



+1 AB 



— (sinh 1 77) — £ i/(sinh x ?7) 
at 



(B13) 



With these expressions at hand, it is a reasonably straightforward exercise to substitute 
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them into (jB4|) . We thereby obtain: 



Jf {(<5P«)(* A fty) - (5aP ij )(5^)} (B14) 
6 { HA ~ I 6 + {X ° VA + XA)3A ] + 8^ Jj X -\ 



- [ {5{\ N)U + 8{\ i + \ V i )'H i } 

+ 8^G J s W&VhiC) - ((5e)(* A lnC)} - fi^T^ 

+ ^ ^ {^eS^ftAB + + 8{aC)0 {n ye + (<5a)e0 (/) + 2<f(are0 (/) ) - 5 (^j J 

+^ / {(5Asmh-S)(5e)-5(sinh-S)(5 A e)} 
87TG y B 

In the above, we have introduced the notation to signify the bulk part of the Hamiltonian, 
namely 

H A = f ^x^oiNn + Vm^ + XUi] . (B15) 

We can now interpret the terms arising in the expression (jB14|) one at a time. The first 
term is an exact variation, which is precisely what is required. This will be part of the 
Hamiltonian which generates the evolution along A. The next term consists of variations 
of the lapse and shift multiplying the bulk constraints. This will clearly vanish provided 
the bulk constraints are satisfied. The remaining terms are all boundary terms. On each 
boundary, we have a contribution containing terms involving the variation of the boundary 
volume forms e and e. From these, we identify the momenta conjugate to the volume forms 
as: 

P = ^-lnC and P = — — sinh" 1 rj . (B16) 
8itG 8ttG 

The equations of motion then guarantee that the evolution of these quantities along A is 
given by £ a as expected. 

Finally we come to the remaining boundary terms at the horizon and outer boundary. 
In all cases, these terms either vanish (because 6^ = 0) or the quantity which appears 
inside variation is held "fixed" by our boundary conditions. Although these quantities are 
held fixed, with the weak boundary conditions introduced in Section QllDl their values will 
change under diffeomorphisms of the boundary. Thus, variations of these quantities do not 
vanish, but are equal to Lie derivatives along a direction X 1 tangent to the boundary surface. 
For a three-dimensional Hamiltonian formulation this means that 



as discussed in f!4.5|) . 
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In addition, we require that the diffeomorphism constraint be satisfied at both the horizon 
and outer boundary. The analysis of the outer boundary was performed in detail in We 
refer the interested reader to that paper (and in particular Appendix C); here we will simply 
use the result to rewrite the terms arising in (|B14|) . Therefore, we turn our attention to the 
horizon. On S, the diffeomorphism constraint can be rewritten in a form (|3.2(Jj) which is 
directly applicable for our current calculation (|B14|) . The constraint equation on the horizon 
allows us to immediately cancel all the terms other than the uo A ones. These can be rewritten 
using: 

x A C v (eu A ) = I S5X A ^-(eu A ) + V A 6 x {eu> A )\ , (B18) 



Additionally, we must express the 5X A term in our preferred co-ordinate system. As in |l5j, 
we can make use of the relation 

5X a = (5X) A (dJ a ) - X A d A {56 a ) , (B19) 

to rewrite 

/ {5X A )e£o A + X {S9 A )C T (euj A ) = [ (5X a )eu a + §6 a £ A {e£o a ) . (B20) 
Js Js 

Finally combining everything from above, we can write the bulk symplectic structure as: 

/ {(SP ij )(S A hij) - (S^XShij)} (B21) 

= s (** - L { + e ^ yA + aA) ^} + L rxA ^ 

- [ (5(A iV)?i + 5(A J + A o r)^} 
Jt, 

+ J { (5 A P) (5e) - (6P) (5 A e)}+J s { (5 A e) (5P) - (<5e) (5 A P) } 
+ [ lJ^ekn + ej A V A \+5t [ C A I J^ek n + ej A V A 



^ 8ttG " JA J J B { 
+ J {(5X a )ej a + Se a 6 A (ej a )} ~^f s {(^ a )^ a + Se a 5 A (eCo 



This is the final desired form. 
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